Unconventional plasmon-phonon coupling in graphene 
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We predict the existence of coupled plasmon-phonon excitations in graphene by using the self- 
consistent linear response formalism. The unique electron-phonon interaction in graphene leads 
to unconventional mixing of plasmon and optical phonon polarizations. We find that longitudinal 
plasmons couple exclusively to transverse optical phonons, whereas graphene's transverse plasmons 
couple only to longitudinal optical phonons. This coupling can serve as a magnifier for exploring 
the electron-phonon interaction in graphene, and it offers novel electronical control over phonon 
frequencies. 

PACS numbers: 73.20.Mf,73.22.Lp,63.22.Rc,78.67.Wj 

The interaction of electrons and crystal lattice vibra- 
tions (phonons) has fundamental implications on prop- 
erties of materials and leads to diverse many-body phe- 
nomena such as superconductivity and charge-density 
waves. The electron-phonon interaction takes an unusual 
form in graphene, a recently discovered two-dimensional 
(2D) material [l[ (see, e.g., (H for a review), and its im- 
plications arc far from being explored. These include 
the breakdown of the Born-Oppenheimer approximation 
Q , anomaly of the optical phonon Q , and nonadiabatic 
Khon anomaly jsj. However, the interaction of collec- 
tive electron excitations (plasmons) and optical phonons 
has not yet been presented for graphene. Plasmons in 
graphene are of fundamental scientific interest 
but they also hold potential for technological applica- 
tions (e.g., in the context of plasmon lasers [9| and meta- 
materials [13 )• Besides the ordinary longitudinal plas- 
mons (transverse magnetic modes) Il2{ , graphene 
also supports unusual transverse plasmons (transverse 
electric modes) @. The hybridization of plasmon and 
phonon modes is a striking manifestation of the break- 
down of the Born-Oppenheimer approximation, because 
it occurs when phonons and electrons arc on compara- 
ble energy scales. Plasmon-phonon coupling has been 
studied in bulk semiconductors Il3 ^ . ll4i , systems with re- 
duced dimensionality (see e.g. [ISMlq ). and in the con- 
text of graphene, plasmons were shown to couple to sur- 
face optical phonons of the substrate (e.g., SiC, which 
is a polar material) [ll|, [l^ [23|. However, to the best 
of our knowledge, the phenomenon of plasmon-phonon 
coupling was not yet presented for an isolated 2D ma- 
terial. Here we predict the coupling of plasmons with 
intrinsic optical phonons in graphene by using the self- 
consistent linear response formalism. We find that, in 
contrast to all other known systems in nature, longitudi- 



nal plasmons (LP) couple only to transverse optical (TO) 
phonons [2l|, while transverse plasmons (TP) couple only 
to longitudinal optical (LO) phonons. The LP-TO cou- 
pling is stronger for larger concentration of carriers, in 
contrast to the TP-LO coupling (which is fairly weak). 
The former could be measured via current experimen- 
tal techniques. Thus, plasmon-phonon resonance could 
serve as a magnifier for exploring the electron-phonon 
interaction, and for novel electronic control (by exter- 
nally applied voltage) over crystal lattice vibrations in 
graphene. 

The low-energy band structure of graphene consists of 
two degenerate Dirac cones at K and K' points of the 
Brillouin zone [Hill [see Fig. [I^a) and (b)], and the 
electron Hamiltonian around K point can be written as 



(1) 



where vp = 10^ ni/s, k = {kx,ky) = —i'V is the wave- 
vector operator, cr = (ax, (Jy), and ax,y are the Pauli spin 
matrices. Wc label the eigenstatcs of Hamiltonian He by 
|s, k) and the appropriate eigenvalues by E^^k — shvpl^-l, 
where s = 1 for the conduction band and s = — 1 for the 
valence band. A technologically interesting property of 
graphene is that the concentration of electrons n, and 
hence the Fermi level Ep = hvpy/irn, can be changed 
via gate voltage [l| . 

The long- wavelength in-plane optical phonon branch in 
graphene consists of two modes (LO and TO) which are 
effectively dispersionlcss and degenerate at energy Huq = 
0A96eV Mm- Let u(R) = [ua(R) - ub(R)]/V2 de- 
note the relative displacements of the sub-lattice atoms 
A and £? of a unit cell specified by a coordinate R [see 
Fig. HJc)]. Then, in the long- wavelength limit R can be 
replaced by a continuous coordinate r and we have 

1 „ 



u(r) - J2 



VNM 



(2) 
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where is the number of unit cells, M is the carbon 
atom mass, q = q{cos 4>q, sin (/)q) is the phonon wave vec- 
tor, fi ~ L,T stands for the polarization, and the po- 
larization unit vectors are e^q = i(cos (/)q, sin(/)q), and 



e^q = i(— sint^q, cos0q). The displacement vector u(r) 
is parallel (perpendicular) to the phonon propagation 
wave vector q for LO (TO, respectively) phonons [see 
Fig. [I{c)]. The phonon Hamiltonian is given by 



ph 



^iq-P/jq 



Qi-ici), 



(3) 



where Q^q and P^q denote phonon coordinate and mo- 
mentum. The electron- phonon interaction takes a pecu- 
liar form in graphene [25| 



-ph 



-V2^a X u(r), 
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(4) 



where cr x u = CTxUy — cr 
carbon atoms distance, 
to write Eq. (g]) as 



b = 0.142 nm is the nearest 
and (3 = 2. Wc find it convenient 



He—ph — L J" ^ ^ jq X e^qQ^q 



(5) 



where jq = —evpL ere 'i"" is the single-particle 
current-density operator, is the area of the system, 

e is charge of the electron, and F = , /^^f . , . 

The electromagnetic field in the plane of graphene 
is completely described by the vector potential A = 
S/jq G^qA^qG*'!'' (scalar potential is gauged to zero, time 
dependence is implicitly assumed, and = L,T denote 
polarizations). The interaction with Dirac electrons is 
obtained by substitution fik — > Kk + eA in Eq. ([T]), which 
leads to 



Hp^r 



evpCT ■ A 



3^q^/^q- 



(6) 



/jq 



By comparing Eqs. (|H) and ^ it follows that electron- 
phonon interaction can be regarded as a presence of an 
effective vector potential 



Acff = F2\{eTqQLci - eLqQTq)e" 



q 



(7) 



k-cff 



that is, He-ph = evpcr ■ Acff. It is evident that Ac 
u(r) = that is the effective vector potential Acff is 
perpendicular to u(r) as illustrated in Figs. [IJc) and (d) 
(see also Rcf. 0), which is responsible for the mixing of 
polarizations in plasmon-phonon coupling. 

As a first pass, let us ignore the phonons and focus on 
the Hamiltonian H = + He-em- Without an exter- 
nal perturbation, the electrons in graphene fill the Fermi 
sea according to the Fermi distribution function fsu- A 
field A^q{uj) oscillating at frequency cj will induce an av- 
erage current density (up to a linear order in the vector 
potential) 

(J^(q,w)) = -Xp(q,w)A^q(w), (8) 

where the current-current response functio n ( including 
2-spin and 2-valley degeneracy) is given by |26| 

/sik /s2k+q 



X^(q,w) =4L' 



SlS2k 



hoj + flWs-,k - hlUs^k+q + irj 

X |(sik|jq-e* |s2k + q)p. (9) 




(b) 
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FIG. 1: (a) Schematic illustration of the lattice structure with 
two sublattices (A and B). (b) The two degenerate Dirac cones 
are centered at K and K' points at the edge of the Brillouin 
zone, (c) A displacement of lattice atoms u(r) is parallel 
(perpendicular) to the propagation wave vector q of a LO 
(TO) phonon. (d) The displacement u(r) creates an effective 
vector potential Aeff perpendicular to u(r) (the sign of Acs 
for the K' point is opposite to that for the K point). 



For the response function y ^(q , uj) we utilize the ana- 
lytical expression from Ref. |27| . The subtlety involved 
with the divergence in Eq. ^ is solved by subtract- 
ing from XL(q,w) [xT(q,'^)] the value Xl{c[,uj = 0) 
[XT(q — > 0,a; = 0)] to take into account that there is 
no current response to the longitudinal [transverse] time 
[time and space] independent vector potential, see [27l[28j 
for details. We would like to note that when working 
with the current-current response function, rather than 
with the density-density response function, the nature of 
the plasmon-phonon interaction (especially the mixing of 
polarizations as shown below) is far more transparent. 

Next, it is straightforward to show from the Maxwell 
equations that an electric current oscillating in a two- 
dimensional plane will induce a vector potential 



(Aiq(cj)) = (JL(q,w)) 



-2a;2eo 



and 



(ATq(w)) = (JT(q,w)) 



Mo 

2v/g2-c^2/, 



(10) 



(11) 



where we have assumed that graphene is suspended in 
air and that there are no other sources present in space. 
This induced vector potential in turn acts on electrons 
in graphene through the interaction Hamiltonian He-em 
which can result in plasmons - self-sustained collective 
oscillations of electrons. From Eqs. ([8|) and (|10|w e get 
the dispersion relation for longitudinal plasmons @, 0, [l2| 



1 



Xi(q,w) = 0. 



(12) 
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From Eqs. ([8]) and pT|) we get the dispersion relation for 
transverse plasmons [8| 



1 



2^q^ - LO-^/c- 



?XT(q,w) = 0. 



(13) 



Longitudinal plasmons are also referred to as transverse 
magnetic modes since they are accompanied by a longi- 
tudinal electric {E) and a transverse magnetic field {B) 
in the plane of graphene. Likewise transverse plasmons 
or transverse electric modes are accompanied by a trans- 
verse electric and a longitudinal magnetic field Q. Dis- 
persion relation of LP (TP) modes is shown by the blue 
dashed line in Fig. [2] (Fig. [3l respectively). Finally 
we note that we arc primarily interested in non-radiative 
modes {q > w/c) in which case fields are localized near 
the graphene plane (z = 0) and decay exponentially: 
E{z),B{z) cx e-l^lV'J'-'^Vc^. 

In order to find the plasmon-phonon coupled excita- 
tions we consider the complete Hamiltonian H = iJg + 
He-em + He-ph + Hph- Wc assumc that the hybrid plas- 
mon phonon mode oscillates at some frequency lo with 
wavevector q (which are to be found). From the equa- 
tion of motion for the phonon amplitudes Q^q one finds 




0.01 0.02 

q /kp 



0.03 



(w^-Wo)(QTq) =i'F(JL(q,w)), 



(14) 



and 



(^2 _ uI){Ql^) = -L^F{JT{q,uj)). (15) 



The electron phonon interaction ([5]) is included as the 
effective vector potential ([7]) in Eq. ([6]), which from Eq. 
^ immediately yields 



FIG. 2: Dispersion lines of hybrid LP-TO plasmon-phonon 
modes (solid lines) and of the uncoupled modes (dashed lines) 
for two values of doping: (a) n = 5 x 10^^ cm^^, and (b) 
71 = 5 X 10^'^ cm~^. The hybridization is stronger for larger 
doping values. Grey areas denote the region of single-particle 
damping. 



The plasmon dispersions relations and ([T5|) appear 
as polos in the Eqs. ([TS]) and ([T^ for the coupled modes, 
which means that the coupling is greatest at the reso- 
nance point where plasmon momentum and energy match 
that of the appropriate phonon mode. We denote this 
point (where the uncoupled plasmon and phonon disper- 
sion cross) by {qc,uJo). One can quantify the strength 
of the coupling effect by calculating the frequency differ- 
ence between the hybrid modes at the wavevector qc in 
units of the uncoupled frequency value: Aw/wq. Finally 
by doping one can change plasmon dispersion which in 
turn changes qc and the strength of the plasmon-phonon 
coupling. 

The dispersion lines for the hybrid LP-TO modes are 
shown in Fig. [2] for two values of doping, (a) n = 
5 X 10^2 cm-^ Ep = 0.261 eV, kp = 3.96 x 10*^ m-\ 
and (b) n = 5 X lO^^ cm^^ Ep = 0.825 cV, kp = 



(JL(q,w)) = XL{q,io){^{ALq{io)) + F{Qtci)), (16) 1-25 x 10^ m \ The strength of the coupling increases 



and 

(JT(q,^)) = XT(q,^)(-(ATq(c^)) - F{Ql^)). (17) 

From Eqs. ()14p - ([T7|) it is clear that transverse (longitu- 
dinal) phonons couple only to longitudinal (transverse) 
plasmons. Apparently, this follows from the fact that LO 
(TO, respectively) phonons are equivalent to oscillations 
of an effective vector potential Ach [see Eq. ([7])], and 
therefore an effective electric field, perpendicular (paral- 
lel, respectively) to q. 

Finally using Eqs. (HOl), and ([16]) we get the 

dispersion relation for the LP-TO coupled mode 



2 2 



^ 2^^eo XL(q,^) 



(18) 



1 



and from Eqs. ([TT]) . p^ . and ([T7)) dispersion relation for 
the TP-LO coupled mode 



L^F^XT(q,c^) 



=XT(q,w)' 



(19) 



with increasing values of doping, and one has for the 
case (a) Aw/wq = 7.5%, and (b) Aw/wq = 15.5%. To 
describe graphene sitting on a substrate (say SiC, which 
is a polar material), one only needs to include the di- 
electric function of the substrate into our calculation. In 
that case plasmons can also couple to surface phonon 
modes of the polar substrate [ill, [l^, [l^l . However, since 
these surface phonons have sufficiently smaller energies 
than optical phonons in graphene out results are qual- 
itatively unchanged in that case. LP-TO hybrid modes 
could be measured by observing the change in the phonon 
dispersion with the Neutron Spectroscopy or Inelastic 
X-ray Scattering. Alternatively, one could use grating 
coupler or Electron Energy Loss Spectroscopy to mea- 
sure the shift in the plasmon energy. Our results imply 
that plasmon-phonon coupling could serve to explore the 
electron-phonon interaction (the frequency shifts at res- 
onance are much larger then those recently measured by 
Raman Spectroscopy Q), and that by externally appling 
voltage one can influence the properties of lattice vibra- 
tions. 

In spite of the fact that the formal derivation of hy- 
brid TP-LO coupled modes is equivalent to the deriva- 
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FIG. 3: Dispersion lines of hybrid TP-LO plasmon-phonon 
modes (solid lines) and of the uncoupled modes (dashed lines) 
for two values of doping: (a) n = 7.5 x 10^^ cm~^, and (b) 
n — 9.5 X 10^^ cm^^. The plasmon-like dispersion is very 
close to the light line q — lj/c; therefore, the ordinate shows 
Aq = q ~ oj/c. 



tion of the LP-TO modes, their properties qualitatively 
differ. First, we note that the dispersion of transverse 
plasmons is extremely close to the light line, and we plot 
Aq ^ q — uj/cvs. frequency w following Ref. For this 
reason, transverse plasmons are expected to have strong 
polariton character and they will be hard to distinguish 
from free photons (also, even a small plasmon lincwidth 
will obscure the distinction). Moreover, they do not exist 
in graphene between two dielectrics with sufficiently dif- 
ferent relative permittivity, where the light lines for the 
dielectrics are separated. Next, transverse plasmons ex- 
ist only in the frequency interval 2Ef > huj > 1.667 Ep 
0, which means that the LO phonon energy must be 
in the same interval for the hybridization to occur. Fig- 
ure |3] shows the dispersion curves of the hybrid TP-LO 
modes for two values of doping, (a) n = 7.5 x 10^^ cm~^, 
Ep ^ 0.101 eV, hp ^ 1.53 x 10* m"!, and (b) n = 
9.5 x 10" cm-2, Ep = 0.114 eV, kp = 1.73 x 10* m^^. 
We observe that the trend here is opposite to that of 
the LP-TO coupling, as the strength of the coupling de- 
creases with increasing doping; specifically, one has for 
the case (a) Aw/wq = 0.17%, and (b) Aw/wq = 0.02%. 
The maximal coupling occurs when 2Ep is just above 
hwQ. and it is zero when hw^ — \.661Ep. We emphasize 
that the strength of the coupling for TP-LO modes is in 
general much weaker than in LP-TO modes. 

Before closing, we note another interesting result which 
is captured by our calculations. Equations (fTS]) and (fT9|) 
for shifts in the energies of TO and LO modes at g = 
reduce to 



1 



(20) 



2u€oc 



which is identical to the result of Rcf. [J], where the 
coupling of optical phonons to single-particle excitations 



was studied, appart from the imaginary term in the de- 
nominator which is zero in [3]. This small but quali- 
tative difference is consequence of phonon coupling to 
the radiative electromagnetic modes, which increases the 
phonon lincwidth. For example, for the doping values of 
?i = 5 X 10^2 cm"^ 5 X 10^^ cm^^^ ^nd 5 x 10^** cm"^ Eq. 
(f20| yields 0.005%, 0.07%, and 0.7%, respectively, for 
the lincwidths, while there is no lincwidth from single- 
particle damping at these doping values. This effect is 
qualitatively unchanged for graphene sitting on a sub- 
strate and could be measured by Raman spectroscopy. 
Finally, we note an interesting solution of Eq. (valid 
for suspended graphene): when the hybrid TP-LO mode 
dispersion crosses the light line it has the same energy 
as the uncoupled phonon mode, i.e., w = wq- In other 
words, LO phonon at a wavevector q = loq/c decouples 
from all (single particle and collective) electron excita- 
tions, while no such effect exists for the TO phonons. 

In conclusion, we have predicted hybridization of plas- 
mons and intrinsic optical phonons in graphene using 
self-consistent linear response theory. To the best of our 
knowledge, this is the first study of such resonance in an 
isolated 2D material. We found that graphene's unique 
electron-phonon interaction leads to unconventional mix- 
ing of plasmon and optical phonon polarizations: longi- 
tudinal plasmons couple exclusively to transverse optical 
phonons, whereas graphene's transverse plasmons couple 
to longitudinal optical phonons; this contrasts plasmon- 
phonon coupling in all previously studied systems. The 
strength of the hybridization increases with doping in LP- 
TO coupled modes, while the trend is opposite for TP- 
LO modes. The LP-TO coupling is much stronger than 
TP-LO coupling, and it could be measured by current ex- 
periments, which would act as a magnifier for exploring 
the clcctron-phonon interaction in graphene. This cou- 
pling is an even more striking example of a breakdown 
of Born-Oppenheimer approximation in graphene than 
the recently measured stiffening of the Raman G peak 
Moreover, plasmon-phonon interaction can serve to 
electronically control the frequencies of lattice vibrations 
in graphene, which could have interesting technological 
implications. In this context we should mention that 
our study opens the way for investigations of plasmon- 
phonon interaction in bilaycr graphene, where phonon 
lasers were recently proposed [29|. 
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